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Abstract 

Scattering of electromagnetic (EM) waves by many small particles 
(bodies), embedded in a homogeneous medium, is studied. Physical 
properties of the particles are described by their boundary impedances. 
The limiting equation is obtained for the effective EM field in the lim- 
iting medium, in the limit a — > 0, where a is the characteristic size of 
a particle and the number M(a) of the particles tends to infinity at a 
suitable rate. The proposed theory allows one to create a medium with 
a desirable spatially inhomogeneous permeability. The main new phys- 
ical result is the explicit analytical formula for the permeability fJ-(x) 
of the limiting medium. While the initial medium has a constant per- 
meability fiQ, the limiting medium, obtained as a result of embedding 
many small particles with prescribed boundary impedances, has a non- 
homogeneous permeability which is expressed analytically in terms of 
the density of the distribution of the small particles and their bound- 
ary impedances. Therefore, a new physical phenomenon is predicted 
theoretically, namely, appearance of a spatially inhomogeneous perme- 
ability as a result of embedding of many small particles whose physical 
properties are described by their boundary impedances. 
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1 Introduction 



In this paper we outline a theory of electromagnetic (EM) wave scatter- 
ing by many small particles (bodies) embedded in a homogeneous medium 
which is described by the constant permittivity eo > 0, permeability fiQ > 
and, possibly, constant conductivity <jq > 0. The small particles are em- 
bedded in a finite domain CI. The medium, created by the embedding of 
the small particles, has new physical properties. In particular, it has a spa- 
tially inhomogeneous magnetic permeability fi(x), which can be controlled 
by the choice of the boundary impedances of the embedded small particles 
and their distribution density. This is a new physical effect, as far as the 
author knows. An analytic formula for the permeability of the new medium 
is derived: 

" (x) = wr 

where 

<Jr(x) = 1 + — ie ujh(x)N(x). 
o 

Here ui is the frequency of the EM field, eo is the constant dielectric pa- 
rameter of the original medium, h(x) is a function describing boundary 
impedances of the small embedded particles, and N(x) > is a function 
describing the distribution of these particles. We assume that in any sub- 
domain A, the number J\f(A) of the embedded particles D m is given by the 
formula: 

M(A) = [ N(x)dx[l + o(l)], a -> 0, 

a J A 

where N(x) > is a continuous function, vanishing outside of the finite 
domain Q in which small particles (bodies) D m are distributed, k £ (0, 1) is 
a number one can choose at will, and the boundary impedances of the small 
particles are defined by the formula 

a K 

where x m is a point inside m— th particle D m , Re h(x) > 0, and h(x) is a 
continuous function vanishing outside f2. The impedance boundary condition 
on the surface S m of the m— th particle D m is E t = Cm[H t , N], where E t 
{H l ) is the tangential component of E (H) on S m , and N is the unit normal 
to S m , pointing out of D m . 

Since one can choose the functions N(x) and h(x), one can create a 
desired magnetic permeability in Q. This is a novel idea, to the author's 
knowledge. 



3 



We also derive an analytic formula for the refraction coefficient of the 
medium in created by the embedding of many small particles. An equation 
for the EM field in the limiting medium is derived. This medium is created 
when the size a of small particles tends to zero while the total number 
M = M(a) of the particles tends to infinity at a suitable rate. 

The refraction coefficient in the limiting medium is spatially inhomoge- 
neous. 

Our theory may be viewed as a " homogenization theory", but it dif- 
fers from the usual homogenization theory (see, e.g., [1] , [2] , and references 
therein) in several respects: we do not assume any periodic structure in the 
distribution of small bodies, our operators are non-selfadjoint, the spectrum 
of these operators is not discrete, etc. Our ideas, methods, and techiques 
are quite different from the usual methods. These ideas are similar to the 
ideas developed in papers [HE], where scalar wave scattering by small bod- 
ies was studied, and in the papers [6], [7]. However, the scattering of EM 
waves brought new technical difficulties which are resolved in this paper. 
The difficulties come from the vectorial nature of the boundary conditions. 
Our arguments are valid for small particles of arbitrary shapes. 

We also give a new numerical method for solving many-body wave- 
scattering problems for small scatterers, see Section 5.2. 

2 EM wave scattering by many small particles 

We assume that many small bodies D m , 1 < m < M, are embedded in a 
homogeneous medium with constant parameters eo, /^o- Let k 2 = uj 2 eo^Q, 
where uj is the frequency. Our arguments remain valid if one assumes that 
the medium has a constant conductivity o"o > 0. In this case eo is replaced 
by eo + i^j- Denote by [E,H] = E x H the cross product of two vectors, 
and by (E, H) = E ■ H the dot product of two vectors. 

Electromagnetic (EM) wave scattering problem consists of finding vec- 
tors E and H satisfying the Maxwell equations: 

V x E = iufi H, V x H = -iue E in D := M 3 \ U A ^ =1 D m , (1) 

the impedance boundary conditions: 

[N, [E, N}\ = Cm[H, N] on S m , 1 < m < M, (2) 

and the radiation conditions: 

E = E + v E , H = H + v H , (3) 
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where C, m is the impedance, N is the unit normal to S m pointing out of 
D m , Eq,Hq are the incident fields satisfying equations ([1]) in all of M 3 . One 
often assumes that the incident wave is a plane wave, i.e., Eq = £e lkax , £ 
is a constant vector, a E S 2 is a unit vector, S 2 is the unit sphere in M 3 , 
a ■ £ = 0, ve and vh satisfy the radiation condition: r(|^ — ikv) = o(l) as 
r := \x\ — > oo. 

By impedance Cm we assume in this paper either a constant, Re £ m > 0, 
or a matrix function 2x2 acting on the tangential to S m vector fields, such 
that 

Re(C m £*,£*)>0 V£*eT m , (4) 

where T m is the set of all tangential to S m continuous vector fields such 
that DivE t = 0, where Div is the surface divergence, and E t is the tan- 
gential component of E. Smallness of D m means that ka <C 1, where 
a = 0.5 maxi< m <M diamZ) m . By the tangential to S m component E l of 
a vector field E the following is understood in this paper: 

E t = E - N(E,N) = [N, [E, N]]. (5) 

This definition differs from the one used often in the literature, namely, from 
the definition E l = [N, E] . Our definition ([5]) corresponds to the geometrical 
meaning of the tangential component of E and, therefore, should be used. 
The impedance boundary condition is written usually as 

£* = C[#*,iV], 

where the impedance £ is a number. If one uses definition ([5]), then this 
condition reduces to ©, because [[N, [H,N]],N] = [H,N\. 

Lemma 1. Problem ([!])-([!]) has at most one solution. 

Lemma [1] is proved in Section 2. 
Let us note that problem ©-([l]) is equivalent to the problems ([6]), ((ZJ), ([3]), 
(|3J), where 

V x V x E = k 2 E in D, H = V X E , (6) 

[N, [E, AH] = x E, N] on S m , 1 < m < M. (7) 

Thus, we have reduced our problem to finding one vector E(x). If E(x) is 
found, then H = Y^. 
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Let us look for E of the form 

E = E +Y j V x / g(x, t)a m {t)dt, g(x, y) = — (8) 

^-^ o 4:ir\x — y\ 

m=l J a ™ 



where t 6 5 m and dt is an element of the area of S m , a m (t) € T m . This E 1 
for any continuous a m (t) solves equation ([6]) in D because Eq solves ([6]) and 

VxVxVx / g(x,t)a m (t)dt = VV • V x / g(x,t)a m (t)dt 

- V 2 V x / g(x,t)a m (t)dt 

J Sm 

= k 2 \7 x / g(x,t)a m (t)dt, x £ D. 

J Sm 

(9) 

Here we have used the known identity divcurlE = 0, valid for any smooth 
vector field E, and the known formula 

- V 2 g(x,y) = k 2 g(x,y) + 5(x-y). (10) 

The integral f s g(x,t)a m (t)dt satisfies the radiation condition. Thus, for- 
mula (|8|) solves problem ©, (|7|), ([3]), 0, if a m {t) are chosen so that bound- 
ary conditions ([7]) are satisfied. 

Define the effective field E e (x) = E™(x) = Ei m \x,a), acting on the 
m— th body D m : 

E e {x) := E(x) - V x f g(x,t)a m (t)dt:= Ei m \x), (11) 

where we assume that x is in a neigborhood of S m , but E e (x) is defined 
for all x € M 3 . Let x m £ Z) m be a point inside D m , and d = d(a) be the 
distance between two neighboring small bodies. We assume that 

lim = 0, lim d(a) = 0. (12) 

We will prove later that E e {x,a) tends to a limit E e (x) as a — > 0, and 
E e (x) is a twice continuously differentiable function. To derive an integral 
equation for a m = a m (t), substitute 



E = E e + Vx J g{x,t)a m (t)dt 

J Sm 
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into (J7|), use the formula 

Cm(*) 



[AT,Vx / g(x,t)a m (t)dt] T = / [N„, [V x g(x,t)\ x=s , a m (t)]]dt ± o 

" Sjn J Sm- 

(13) 

(see, e.g., [3]), the -(+) signs denote the limiting values of the left-hand 
side of (|13p as x — )• s from L> (D m ), and get the following equation (see 
Appendix): 

a m (t) = A m a m + f m , 1 < m < M. (14) 

Here A m is a linear Fredholm-type integral operator, and f m is a contin- 
uously differentiable function. Let us specify A m and f m . One has (see 
Appendix): 

f m = 2[f e (s),N s ], f e (s):=[N s ,[E e (s),N s ]]-^-[V xE e ,N s }. (15) 
Condition ([7]) and formula (fl~3j) yield 

/e(s) + + [ / Ws, [V s g(s,t),a m (t)]]dt,N s 



W -[VxVx / g(x,t)a m (t)dt,N s }\ x ^ s = 



(16) 



(17) 



^Mo Js„ 
Using the formula V x Vx = graddiv — V 2 , the relation 

V^Vz • / g(x,t)a m (t)dt = V X (-V t g(x,t),a m (t))dt 

= V X J g(x,t)Biva m (t)dt = 0, 
where Div is the surface divergence, and the formula 

-V 2 / g(x,t)a m (t)dt = k 2 J g(x,t)a m (t)dt, x G D, (18) 

J Sm " Sm 

where equation (jlOp was used, one gets from (|16p the following equation 

- [N s ,a m (s)] + 2/ e (s) + 2Ba m = 0. (19) 

Here 



Ba m :=[l [N s ,[V s g(s,t),a m (t)]]dt,N s ] + ( m iue [l g(s,t)a m (t)dt,N s }. 

(20) 



Take cross product of iV s with the left-hand side of (|19p and use the formulas 
N s ■ a m (s) = 0, f m := f m (s) := 2[f e (s),N s ], and 

[N s ,[N s ,a m {s)}} = -a m {s), (21) 

to get from ([H?]) equation (TH|) : 

a m (s) = 2[f e {s),N s ] - 2[N s ,Ba m ] := A m a m + f m , (22) 

where A m a m = —2[N S , Ba m ]. The operator A m is linear and compact in 
the space C(S m ), so that equation ([22]) is of Fredholm type. Therefore, 
equation ()22[) is solvable for any / m € T m if the homogeneous version of 
(|22p has only the trivial solution <r m = 0. In this case the solution a m to 
equation (j22l) is of the order of the right-hand side f m , that is, 0{a~ K ) as 
a — ?> 0, see formula (|15p . Moreover, it follows from equation (j22|) that the 
main term of the asymptotics of a m as a — > does not depend on s € £ m . 

Lemma 2. Assume that a m G T m , <r m £ C(S m ), and a m (s) = A m a m . Then 
a m = 0. 

Lemma [2] is proved in Section 2. 

Let us assume that in any subdomain A, the number A/"(A) of the em- 
bedded bodies D m is given by the formula: 

AT(A) = / N(x)dx[l + o(l)l, a^O, (23) 

a 2 K Ja 

where iV(a;) > is a continuous function, vanishing outside of a finite domain 
in which small bodies D m are distributed, n € (0, 1) is a number one can 
choose at will. We also assume that 

_ h{x m ) 

where Re h(x) > 0, and /i(x) is a continuous function vanishing outside 
Let us write ([8]) as 

M Af . 

E(x) = E (x)+^2[V x g(x,x m ),Q m }+^2 Vx / (g(x,t)-g(x,x m ))a m (t)dt, 

771=1 777=1 , ' 5 ' m 

(25) 

where 

Q m := [ a m {t)dt. (26) 



S 



Since a m = 0(a K ), one has Q m = 0(a 2 K ). We want to prove that the 
second sum in (|25[) is negligible compared with the first sum. One has 

ii := \[V x g(x,x m ),Q m }\ < O (maxj j^}) °( a2 ~ K )> ( 27 ) 

32 ■= |V x J (g(x,t) - g{x,x m ))a m {t)dt\ < aO (max j , ^ | ^ 0(a 2 " K ), 

(28) 

and 

^ = O (max{^,/ca}) ->• 0, ^ = o(l), a ->• 0. (29) 
Thus, one may neglect the second sum in (|25p . and write 

M 

E{x) = E (x) + \y*9(?, x m ),Q m ] (30) 

m=l 

with an error that tends to zero as a — > 0. 

Let us estimate Q m asymptotically, as a — > 0. Integrate equation ([22]) 
over 5 m to get 

Qm = 2[ [f e (s),N s ]ds-2 I [N s ,Ba m ]ds. (31) 

We will show in the Appendix that the second term in the right-hand side 
of the above equation is equal to — Q m plus terms negligible compared with 
the first one as a — >• 0. Thus, 



Qm= f [f e (s),N s ]ds. 

J Sm. 



Let us estimate the first term. It follows from equation (|15p that 

[N s , f e ] = [N s , E e ] - S^[N S , [V x E e , N s }]. (32) 

If E e tends to a finite limit as a — > 0, then formula (|32|) implies that 

[N a , fe] = O(Cm) = O f^j , a -> 0. (33) 

By Lemma[2]the operator (7 — ^Ln) -1 is bounded, so cx m = O (^«), and 

Q m = O (a 2 - K ) , a -»• 0, (34) 
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because integration over S m adds factor 0(a 2 ). As a — > 0, the sum ([30 
converges to the integral 



(35) 



E = E Q + Vx [ g(x,y)N(y)Q(y)dy, 
Jn 

where Q(y) is the function such that 

Qm = Q(x m )a 2 - K . (36) 
The function Q(y) can be expressed in terms of E: 

Q(V) = -y%)^eo(V x E)(y), (37) 

see Appendix. Here the factor f appears if D m are balls. Otherwise a 
tensorial factor c m , depending on the shape of S m , should be used in place 
off. 



Thus, equation (]35|) takes the form 
E(x) = E (x) - — iueoV x / g(x,y)V x E{y)h(y)N(y)dy. (38) 



Let us derive physical conclusions from equation (|38|) . Taking V x Vx of 

{BSD y ields 



VxVx£ = /c 2 ^^) 

- ^iwe V x (grad div - V 2 ) / g{x,y)V x E(y)h(y)N(y)dy 
6 Jn 

= k 2 E - k 2 ^-ioje V x / g(x,y)X7 x E(y)h(y)N{y)dy 
3 Jq 



— iwe V x (V x E(x)h(x)N(x)) 
3 

fc 2 £;(x) iuje h(x)N(x)V x V x E 

3 

-f iw€ [V(/»(a;)JV(a;)), V x E(x)\. 



(39) 



Here we have used the known formula V x grad = 0, the known equation 
(|10p . and assumed for simplicity that h(x) is a scalar function. It follows 
from (EH) that 



§T Ll 

VxVxB = K 2 (x)E 3 "7 [V(h(x)N(x)),VxE(x)}, (40) 



3 -iueo 
1 + tfiue h(x)N(x) 
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where 

k 



K\x)= ~ k 2 = u 2 e n . (41) 

1 + ^-iujeoh{x)N (x) 



UVxE = iufx(x)H and V x H = —ioje(x)E, then 



V x V x E = u 2 e(x)u(x)E + [ V/ ^ V x E}. (42) 

u(x) 



Comparing this equation with (|40p , one can identify the last term in (|40p as 
coming from a variable permeability fi(x). This /u(x) appears in the limiting 
medium due to the boundary currents on the surfaces S m , 1 < m < M. 
These currents appear because of the impedance boundary conditions ([T]). 
Let us identify the permeability //(x). Denote \&(x) := 1 + ^-icoeoh(x)N(x). 
Let e(x) = 6q, eo = const, and define /i(x) := -^y- Then ivT 2 = w 2 eo/x(a;), 

and = — . Consequently, formula (jlOj) has a clear physical mean- 

ing: the electromagnetic properties of the limiting medium are described by 
the variable permeability: 

u ( x ) = = /fo m 



3 Conclusions 



The limiting medium is described by the new refraction coefficient K 2 {x) 
( see (|4ip ) and the new term in the equation (140 j) . This term is due to the 
spatially inhomogeneous permeability n(x) = generated in the limiting 
medium by the boundary impedances. The field E{x) in the limiting medium 
( and in equation (j4"0|) ) solves equation ([38]) . 

Therefore, we predict theoretically the new physical phenomenon: by em- 
bedding many small particles with suitable boundary impedances into a given 
homogeneous medium, one can create a medium with a desired spatially in- 
homogeneous permeability (|43p . 

One can create material with a desired permeability fj,{x) by embedding 
small particles with suitably chosen boundary impedances. Indeed, by for- 
mula (I43p one can choose a complex-valued, in general, function h(x), and 
a non-negative function N(x) > 0, describing the density distribution of 
the small particles, so that the right-hand side of formula (l43l) will yield a 
desired function (a(x). 
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4 Proofs of Lemmas 

Proof of Lemma \Ji 

From equations (pQ) one derives (the bar stands for complex conjugate): 

/ (H -V x E - E -V xH)dx = (iufj, \H\ 2 - iuje \E\ 2 )dx, 
Jd r Jd r 

where Dr := DHBr, and R > is so large that D m C Br := {x : \x\ < R} 
for all m. Recall that V • [E, H] = H-X7xE-E-VxH. Applying the 
divergence theorem, using the radiation condition on the sphere Sr = 8Br, 
and taking real part, one gets 

M . M . 

= ^ Re / [E,H] ■ Nds = Re / C^T^ • E^ds, 

m —l J S m m —l J S m 

where Ef is the limiting value of E t on S m from D, E = £ m [H,N]. This 
relation and assumption imply E^ = on S m for all m. Thus, E = H = 
in D. 

Lemma [1] is proved. □ 
Proof of Lemma [H 

If o~ m = A m a m , then the functions H = and E(x) = Vx f s g(x, t)a(t)dt 
solve equation ([1]) in D, E and H satisfy the radiation condition, and , con- 
dition ([2]). Thus, E = H = in D. Consequently, 

= V x V x / g(x,t)a m (t)dt = (grad div - V 2 ) / g(x,t)a m (t)dt 

= k 2 g(x,t)a m (t)dt, x G D. 
This implies cr m (s) = 0. 

Lemma [2] is proved. □ 



5 Appendix 

In Section 5.1 equation (|38p is derived. In Section 5.2 a linear algebraic 
system (LAS) is derived for finding vectors Q m in equation (j36|) . 
5.1. Boundary condition ([7]) yields 

= [N[E e ,N]} - ^[V x E e , N] + [N, [V x f g(s,t)a m (t)dt, N}} 
'-7'u Vs... 

_J^[VxVx f g(x,s)a m {t)dt,N}. 



12 



Let us denote 

f e := [N, [E e ,N]} -^[Vx E e ,N}. 
One has V x Vx = curlcurl = graddiv — A, and 

V x -/ g[x,t)a m (t)dt = - I (V t g(x,t),a m (t)) dt = g(x,t)Vfcr m (t)dt = 0, 

J S m. J S m. J S m. 



and 



-V x / g(x,t)a m (t)dt = k / g(x,t)a m (t)dt, 



because — V^g(x,t) = k 2 g(x,t), x ^ t, see (fTUj) . Thus, using (fT3j) . one gets: 
= f e + [f [N s , [V s g{s,t),a m (t)}}dt,N s ] + ^[a m (s),N s ] 

J Srn, 

+ — [iV s , / ^(s,t)o- m (t)cft]. 



Cross multiply this by N s from the left and use the relation N s ■ a m {s) = 0, 
to obtain 

= [N s , f e ] + [N„ [ I [N„ [V s g(s, t), a m (t)]]dt, N s }] + ^a m (s) 



S;-i 



- ( m iueo[N s , [N s , I g(s,t)a m (t)dt}}. 
Note that 



[N„[ [N s ,[V s g(s,t),a m (t)}]dt,N s }}= [N s ,[V s g(s,t),a m {t)]]dt 

J Sm " Sm 



N s I ([N s ,[V s g{s,t),a m {t)}}dt,N s ), 

Sm 

[N s ,[V s g(s,t),a m (t)])dt. 



where the integral before the last one vanishes because its integrand vanishes 
as the dot product of two orthogonal vectors. Consequently, 



<r m (t) = 2[f e (s),N s ] + 2( m iue [N s , [N s , / g(s,t)a m (t)dt}} 

J S m 

-2 [N s ,[V s g(s,t),a m (t)]]dt := Aa m + f m , 

J Sm 
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which is equation (fT4"|) . and 

f m :=2[f e (s),N s ], 



which is equation (fT5"]) . 
Denote 

Qr, 



/ a m (s)ds. 

J Sm. 



One has 



[[N s ,[E e (s),N s ]],N s }ds = f [E e {s),N s ]ds 



V ' x x E e dx, 



and 



/ [[V xE e ,N s ],N s ]ds = -( f VxE e ds- f N S (V x E e , N s )ds) 

f Aira 2 
= - V x E e ds + — — V x B e (x ra ) 

/ C O 



87ra 2 „ , . 
= ^— V x E e (x m ), a^O. 

Here we have used the formulas 



/ V x E e ds ~ 47ra 2 V x E e (x m ), a ->• 0, 

and 

47rct 2 

NiNjds = —z—6ij, 
s 6 

where 5 is a sphere of radius a, {iVi}? =1 are Cartesian components of the 
outer unit normal to the sphere S, and 5ij = if i ^ j, 8a = 1. 
Thus, 

f m (s)ds r ^-( m a 2 iuje V x £ e (x m ) = 0(a 2_K ), 



as a — > 0, provided that 



Cm = %l, 0<K<1. 



Let us now show that the term f e Aa m ds contributes — Q, so 

Q m ~ 0.5 / f m (s)ds, a^O. (44) 
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One has 



2 ds [N s ,[V s g(s,t),a m (t)]]dt 

" Sm " Sm 

-2 I ds j dt (v s g(s,t)(N s ,a m {t)) - a m (t) ) dl 



Since 



-2 j ds j dtV s g(s,t)(N s ,a m {t))+ f a m {t)dt2 j ds d ' J[ " h 

J Sm J Sm " Sm " Sm 

J S m s J Dm 



one gets 

I := j dta m (t)2 f ds 99 ^ = - f a m {t)dt-2k 2 f dta m {t) f dxg{a 

J Sm " Sm & ** $m JSm JD 

rn 

Therefore 

I ■= —Qm + h, 

where the term I\ is negligible compared with I. 
If J s \(Tm(t)\dt < oo and j s a m (t)dt ^ 0, then 

a m (t)dt\ > [ / dta m (t) I dxg{x,t)\, 

Sm J Sm " Dm 

because | J D dxg(x,t)\ = 0(a 2 ) if x G D m . 
One has 

| - 2 / ds dtV s g(s,t)(N s ,a m (t))\ < | / a m (t)dt\ = \Q m \, 

J Sm " Sm " Sm 

because \(N s ,a m (t))\ = 0(\s — t\) as \s — t\ — > 0. Therefore, 

Q m = 0.5 / a m (t)dt — C m a 2 iue Q V x E e (x m )ds, a -> 0. (45) 

JSm 6 

This yields the following formula (cf (|30p ): 

8vr M 

E(x) = Eo(x)-—iue ^( m a 2 [Vg(x,x m ),VxE e (x m )], a -> 0. (46) 
■j , 

m=l 
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or, 

M 

r . 



87T 

£(x) = £ o (a0 - — iue ^ /i(x m )a 2 ~ K [V x g(x, x m ), V x B e (xm)] • (47) 

m=l 



Passing to the limit a — > 0, one obtains 

E e = Eq(x) - —ioje / [X7 x g(x, y), V x E e (y))h(y)N(y)dy, (48) 

where /i(x) is the function in the formula £ m = , and iV(x) is the 

function in the definition of 7V(A). The above passage to the limit is done 
by Theorem 1 from [7J, p. 206. It uses the convergence of the collocation 
method for solving equation ()38f) . see [6]. Writing E e = E for the limiting 
field yields equation (j38|) , 

5.2. In this Section a numerical method is developed for solving many- 
body wave scattering problem when the scatterers are small in comparison 
with the wavelength. The method consists of a derivation of a linear alge- 
braic system for finding vectors V m := (V x E)(x m ), 1 < m < M. If V m are 
found, then by formulas ([37]) and ([36]) one finds 



Qm = — —mu}e h(x m )a 2 ~ K V m , (49) 



and, by formula (j30|) . the field E(x). 

Let us derive linear algebraic system for finding V. t 



in ■ 



Apply Vx to equation ([30|) . let x = Xj, 1 < j < M, and replace Yl m =\ 
by the sum £m^,m=i- 
Then one obtains 

8vr M 

='Poj--^-iue a 2 ~ K ^ (graddw-V 2 )g(x,x m )\ x=X:j h(x m )V m , (50) 

7717^', m=l 

where 1 < j < M, and 

Poi := (V x E )( Xj ). (51) 

Equation (|50p is a linear algebraic system for finding V m . In the above 
derivation we have used the formula 

V x [A, B] = (B, V)A - (A, V)B + A(V, B) - S(V, A), 

which yields 

V x [V 5 , Q m ) = (Qm, V)V 5 - Q m V 2 <?, 
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where 



33 ^2, 



(Qm,V)Vg = ^2 e i^2Qmi' Q x g xi , Qmi = (Qm,ei), 



1=1 i'=l 



and ej, i = 1,2, 3, is the orthonormal Cartesian basis of 
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